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We present a compelling evidence for the opening of a bandgap in exfoliated bottom-gated bi- 
layer graphene by fitting the gate-voltage modulated infrared reflectivity spectra in a large range of 
doping levels with a tight-binding model and the Kubo formula. A close quantitative agreement be- 
tween the experimental and calculated spectra is achieved, allowing us to determine self-consistently 
the full set of Slonczewski-Weiss-McClure tight-binding parameters together with the gate-voltage 
dependent bandgap. The doping dependence of the bandgap shows a good agreement with the 
existing calculations that take the effects of self-screening into account. We also identify certain 
mismatches between the tight-binding model and the data, which can be related to electron-electron 
and electron-phonon interactions. 



I. INTRODUCTION 

Bilayer graphene has recently attracted much atten- 
tion motivated by a broad spectrum of unusual elec- 
tronic properties and a number of possibilities for ap- 
plications. It represents the simplest system, where 
the effects caused by a coupling between graphene lay- 
ers can be studied and exploited. Although the inter- 
layer coupling is much weaker than the in-plane chem- 
ical bonding, it results in profound differences between 
electronic and transport properties of monolayer and bi- 
layer graphene , as exemplified by the anomalous quan- 
tum Hall effectP^m. Another notable dissimilarity is re- 
lated to the behavior of these systems in a perpendicu- 
lar electric field. While in zero field both of them are 
zero-gap semiconductors (or zero-overlap semimetals), in 
bilayer graphene a bandgap is generated in the presence 
of the field, due to an introduced as ymmet ry of the elec- 
trostatic potential on the two planed^'^'^. Importantly, 
the bandgap can be tuned continuously, either by apply- 
ing gate voltage* or chemically^, which, in combination 
with a high mobility of charge carriers, opens new unex- 
plored avenues for using bilayer graphene in field effect 
transistors (FETs)ff2IllI12l and other electronic devices. 

Using angle-resolved photoemission spectroscopy 
(ARPES), Ohta et al.l^ indeed observed a band gap m 
bilayer graphene epitaxially grown on top of silicon 
carbide and doped chemically with potassium. A more 
" clean" way of introducing charge carriers is by applying 
electric field using gate electrodes. This technique can 
be most easily applied to exfoliated samples, produced 
by micromechanical cleavage of graphite^"^-. Apart from 
the bandgap generation, applying a gate voltage has 
also the usual doping effect. In order to control the 
doping and the bandgap independently, Oostinga et 
fabricated two electrodes on both sides of the sample and 
found an insulating state, when gate voltages of opposite 



sign were applied to the electrodes. This showed the 
existence of the bandgap, although the determination 
of its exact value from the DC measurements was not 
possible. 

Infrared spectroscopy, which is one of the most 
direct methods to measure the bandgap and other 
band characteristics in conventional semiconductors, is 
clearly technique of choice also in the case of bilayer 
graphene. The two-dimensionality of this material per- 
fectly matches geometrical requirements of an optical 
experiment. Moreover, the possibility of changing the 
chemical potential and the bandgap with the gate volt- 
age suppli es an unprecedented amount of additional 
informatiorj^ i-'^ ' i6 | i7 | i8|i9 | 20 | 2 il22l compared to standard 

optical measurements. 

The tight-binding theory is widely used to describe 
the low energy tt bands in graphitic materials. In the 
case of graphite, a set of tight-binding parameters, known 
as the Slonczewski-Weiss-McClure (SWMcC) modePl2ll^ 
was very successful in describing quantitatively the de 
Haas-van Alphen effect and optical spectreP^. There- 
fore we one can expect that it will also apply to bilayer 
graphene, if the bandgap is properly included. It appears 
that all SWMcC parameters influence the optical con- 
ductivity for photon energies below 1 eV. However, the 
effects of different parameters are rather dissimilar and 
not all of them can be easily extracted from the spectra. 

Several calculations of the optical conductivity of bi- 
layer graphene within the tight binding approach were 
done. Nilsson et al.^^ and Abergel and FalkcP3 consid- 
ered the simplest model, which contains only the nearest- 
neighbor in-plane and interplane hopping terms (70 '^-^ 3 
eV and 71 ~ 0.4 eV respectively) and found that the opti- 
cal conductivity is marked by a profound structure at the 
photon energy huj ~ 71 (we shall refer to this structure as 
the 71 - peak). Nicol and Carbott^^H extended this model 
to include a bandgap and finite doping. Zhang et dX^^ 
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studied the role of additional parameters responsible for 
the electron-hole asymmetry. 

Using infrared techniques, Wang et aP^ observed, in 
agreement with aforementioned predictions, a profound 
anomaly in the spectra of gated bilayer graphene at 
/ttj ~ 7i. Later on, a clear electron-hole asymmetry was 
found in the infrared spectra by Li et aP^ and by Kuz- 
menko et aP^, which allowed determining more SWMcC 
parameters. However, the presence of the bandgap in 
these studies, although expected, was not evident. In 
Ref ll9l no experimental signatures of the bandgap were 
reported and in Ref 1201 onlv a partial agreement between 
the experimental data and a tight-binding calculation in- 
cluding a bandgap obtained theoretically^ was found. As 
it will be discussed below, this is explained by the fact 
that the manifestation of the bandgap in doped graphene 
is more subtle, and therefore requiring more accurate op- 
tical measurements and delicate analysis than in the case 
of an undoped sample. In the latter case a sharp absorp- 
tion threshold corresponding to the electron-hole excita- 
tions across the gap is expected. 

Such a structure was indeed observed by Zhang et al.l^, 
who measured infrared absorption of double gated bilayer 
graphene, where the electric field and doping could be 
controlled independentljEI. A bandgap up to 250 meV 
was observed in an undoped sample in the presence of 
the largest applied fields, which is a rather promising 
sign for the use of bilayer graphene in electronics. Al- 
though for the large values of the bandgap the match 
between the experiment and a tight-binding calculation 
was very close, at low gate voltages, where the absorp- 
tion threshold was beyond the experimentally accessible 
range, the quantitative agreement between the experi- 
mental and theoretical curves turned out to be rather 
poor. 

In general, the quantitative agreement between in- 
frared spectra of bilayer graphene and tight-binding 
model was up to now not very good. Apart from the men- 
tioned discrepancies, the measured height of the 71 peak 
in RefflSl was about two times larger than the calculated 
value. Mak et al. reported the opening of a bandgap 
in top-gated bilayer graphend^, where a small thickness 
of the gate insulating layer allowed very efficient doping. 
The quantitative agreement with the tight-binding the- 
ory, however, was limited, which the authors attributed 
to many-body correlation effects. 

A common problem that one encounters when analyz- 
ing infrared spectra of graphene is their sensitivity to 
several band parameters, including the bandgap, making 
their separate extraction quite complicated. Another is- 
sue is a possible inhomogeneity of the doping level, which 
has a similar effect on the optical spectra as an elevated 
temperature, as will be shown below. In the double- 
gate experiments, the deposition of the top gate on top 
of graphene may affect the band structure and increase 
electronic scattering, not to mention a more complicated 
optical multilayer model that has to be used in order to 
extract the optical conductivity of graphene. Therefore, 



in the present study we fit directly the measured refiectiv- 
ity spectra of bottom gated bilayer graphene with a tight- 
binding model that involves the bandgap, the SWMcC 
parameters, scattering rate, temperature and the im- 
purity doping as adjustable parameters. We find that 
a good quantitative agreement can in fact be achieved, 
which implies that the band structure of bilayer graphene 
is well captured by the tight binding model. Neverthe- 
less, certain discrepancies remain that may eventually be 
related to many-body effects. 

The remaining sections are organized as follows. In 
section |ll] we describe the sample preparation, infrared 
experiment and a specially developed technique of direct 
fitting of the whole set of reflectivity spectra with a tight- 



binding model and the Kubo formula. In section III the 
results of the infrared measurements, their fits, and the 
doping dependence of the bandgap and the chemical po- 



tential are shown. In section IV we discuss the reliability 
of extracting the bandgap and compare its gate-voltage 
dependence with the existing theoretical predictions. We 
demonstrate a practical way of visualizing important fea- 
tures of the band structures using the experimental re- 
fiectivity data. We also discuss the advantages and limi- 
tations of the tight-binding model in describing infrared 
spectra. 



II. TECHNIQUES 
A. Sample 

A relatively large (~ 100 /im) fiake of bilayer graphene 
produced by micromechanical cleavag^I^l of graphite sin- 
gle crystals on top of a Si02 (300 nm)/n-Si substrate 
was chosen for infrared experiments. Lithographically 
deposited leads to the fiake and the gate (doped silicon) 
allowed a simultaneous measurement of the DC conduc- 
tivity ctdc and the infrared refiectivity R as functions of 
the gate voltage Vg. By sweeping Vg from the positive 
to the negative values, one continuously varies the dop- 
ing from electron to hole type. The charge concentration 
can be determined using the relation n — a(Vg — Vcn)i 
where the coefficient a — 7.2 x 10^° cm~^/V is given by 
the electric capacitance of the oxide layer. Usually the 
charge neutrality point Vcn is related to the minimum 
Vmin of the (JDciVg) curvc. However, from the fits of 
the optical spectra described below we found that Vqn 
(which is about -22 V for the presented series of mea- 
surements) is slightly different from Vmin (~ -29 V). It 
is likely that this is related to the large flake dimensions 
and a doping inhomogeneity near electrical contacts. 

The negative value of Vcn is due to charge transfer 
by adsorbed gas molecules (environmental doping) In 
order to reduce this effect, the sample was annealed in a 
H2-N2 atmosphere at 150 °C before each series of mea- 
surements. Although the annealing indeed shifts Vcn 
closer to zero bias, we found that, in contrast to the case 
of monolayer graphene, it never results in Vcn = 0. It is 
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possible that the remaining dopants are located either be- 
tween the flake and the substrate or are even intercalated 
between the carbon layers. As we will see below, this cor- 
relates with the doping dependence of the bandgap. 



B. Optical experiment 
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FIG. 1: (Color online) (a) Absolute (gold-normalized) reflec- 
tivity of bare oxide and graphene at Vg = -1-100, -20 and -100 
V). The substrate temperature is 10 K. (b) Reflectivity of 
graphene normalized to bare oxide, (c) Self-normalized re- 
flectivity of graphene Rgr{Vg) / Rgr{~20 V). 



mounted in a He flow cryostat; a specially made sample 
holder allowed us inserting the sample with a minor delay 
after the annealing. 

Absolute reflectivity spectra of the bare substrate 
Rsubi^) and the flake Rgr{uj) were obtained using a gold 
patch deposited near the sample as a reference (FigjTjL). 
The substrate and graphene spectra look similar: they all 
show a prominent minimum at 0.75 eV due to the Fabry- 
Perot effect in the Si02 layer and intense peaks below 
0.15 eV originating from dipole-active lattice vibrations 
in silicon oxide. The change of reflectivity introduced 
by graphene is better seen in FigjljD, where Rgr{u;) is 
normalized to Rsubioj). One can see that graphene intro- 
duces a significant infrared contrast, especially close to 
the Fabry-Perot minimum. Notably, the same effect in 
the visible range makes graphene detectable by human 
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From FigjTJ) it is also clear that varying the gate volt- 
age has a strong effect on reflectivity. To see it even bet- 
ter, one can normalize Rgr by its value at Vg — —20 V, 
which is close to the charge neural point (Fig|lJ:). Note 
that the Si02 phonon features do not fully cancel after 
the normalization, because of a non-linear character of 
the contribution of the substrate to the spectra. At 0.2 
eV one can see a sharp structure related to the infrared 
active phonon mode in graphene, which has a Fano shape 
and strongly increases as a function of the gate voltage 
due to a coupling to electronic interband transitions, as 
discussed in detail in Ref. |3T] As we shall see below, 
the other changes are due to a combination of the doping 
effect and the opening of the bandgap. 

One should note that the diffraction of electromagnetic 
radiation may affect the measured reflectivity at low en- 
ergies, where the wavelength becomes comparable to the 
spot size. Due to the diffraction and other systematic un- 
certainties, the absolute accuracy of Rgr and Rgr/Rsub is 
about 0.01-0.02 as we determined by varying the position 
and the size of the beam spot on the sample. In contrast, 
the self-normalized reflectivity Rgr{Vg) / Rgr{— 20 V) can 
be measured much more accurately (with the uncertainty 
less than 0.002) since it does not require any reference 
measurement and therefore does not involve any mechan- 
ical movements. In order to minimize the influence of 
weak drifts of the signal, taking spectrum at each gate 
voltage was immediately followed by a separate measure- 
ment at the charge neutral point. 



Optical reflectivity spectra in the range of photon ener- 
gies 0.06 - 1 eV were collected at nearly normal incidence 
using an IR microscope Bruker Hyperion 2000 attached 
to a Fourier transform spectrometer with a standard glo- 
bar source. Appreciable signal could be obtained using 
the beam spot down to 10-15 /im. However, prior to each 
series of measurements we optimized in a try-and-error 
fashion the spot size in order to maximize the signal-to- 
noise ratio while avoiding unphysical spectral artefacts 
related to the finite sample dimensions and apparatus is- 
sues. The spectral resolution was 1 meV. The sample was 



C. Data modelling 

For the data analysis we chose the SWMcCl22i24| tight 
binding description of the tt bands in bilayer graphene^. 
We begin with considering it simply as a band structure 
parametrization in order to extract information about 
the electronic bands from the infrared data. Later on 
we shall discuss the limitations of this description as a 
physical model based on the quality of the obtained fits. 

The structure of the Bernal-stacked bilayer graphene is 
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FIG. 2: (Color online) Crystal structure of Bernal stacked 
bilayer graphene and the considered tight-binding parameters. 



shown in Fig|2] Each layer has two sublattices: Al and 
Bl (bottom layer) and A2 and B2 (top layer). Atoms 
Al and A2 are on top of each other, while the atoms Bl 
and B2 are shifted horizontally by the vectors i5i,2,3 con- 
necting nearest neighbors within one layer. The SWMcC 
Hamiltonian involves the in-layer nearest-neighbor hop- 
ping 7o and three inter layer hopping terms: 71 (between 
Al and A2), 73 (between Bl and B2) and 74 (between 
Al and B2 or between Bl and A2). In addition, the on- 
site energy difference A between the positions Al and 
Bl (A2 a nd B2) i s introduced. Following the standard 
procedur^^'S'^ElM]^ ^^^^ ^^i^ SWMcC model an ex- 
tra parameter U in order to describe the difference be- 
tween the (screened) electrostatic potential of the top 
and the bottom layers in the external field. Note that 
U gives exactly the separation between the electron and 
hole bands at the K point and is slightly larger than the 
true bandgap Ag « \U\^i/{U'^ + 7i)^^^- In the basis 
|B1), |Al), |A2) and |B2) the present Hamiltonian reads 
as follows: 




700 -740 730* 

A 71 -740 

71 A + U 700 

-740* 7o0* U 



(1) 



where 
tum. 



, q is the electronic momen- 



Within the linear response theory, the Kubo formula 
can be used to calculate the complex optical conduc- 
tivity cr(w). In the case of a thin layer, a physically 
more relevant quantity is the optical sheet conductance 
G{u!) = a{uj)d where d is the layer thickness. The total 
conductance consists of the Drude, the interband and the 
high-frequency terms: 



G{co) = Go(c^) + Gib(c^) + Goo(w). 



(2) 
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where Go = e'^/'ih w 6.08 x 10"^ Vt-'^ is the uni 



versal AC conductance of monolayer graphen^^^I^IIMEZl 

and graphite'^*, e^^i are the electronic band energies, 
r is the electronic broadening parameter and /(e) — 

|l-|-exp|^| is the Fermi-Dirac distribution. For 

the derivation of these formulas we can refer, for exam- 
ple, to Ref. 36, where it was done in the case of monolayer 
graphene. 

The term Goo(w), which is purely imaginary in the 
considered spectral range, absorbs contributions from 
all high-frequency core-level and valence band electronic 
transitions, in particular, the ones involving a bands. We 
assume that it does not contribute to the doping depen- 
dence of the optical spectra. 

The chemical potential fi is determined implicitly by 
the doping level via equation: 



1 

2^ 



= n = aiVg - Vcn) (5) 



The first two terms can be obtained using the expressions: 



We subtract 1 /2 in the brackets because the doping level 
is counted relative to half filling of the tt bands. The 
spin degeneracy is already included in equations ([s]), ^ 
and ([5| but the valley degeneracy is not. It appears that 
for the considered range of energies and temperatures it 
is sufficient to perform the momentum integration only 
in a circle of about 1 percent of the total 2D Brillouin 
zone around the K point (and multiply by 2 to account 
for the valley degeneracy). In practical implementation, 
we replace the integration with a summation over ~ 10^ 
q-points. 

In equations ^ and ^ we introduced two different 
scattering rates for the Drude and the interband compo- 
nents. Above 0.1 eV, the real part of the Drude conduc- 
tance is much smaller than the imaginary part, and the 
latter is only very weakly affected by the Drude scatter- 
ing in this range. Since the data is not sensitive to the 
Drude scattering rate, we adopted Fq = 5 meV, which 
corresponds to the value found in graphitd^. We assume 
that the interband scattering F is constant, i.e it is en- 
ergy, momentum and band independent. By doing this, 
we neglected the energy dependent electron-phonon and 
electron-electron scattering processes, which is perhaps 
the most serious limitation of the present model. We will 
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see that it is likely in the origin of some deviations of the 
model curves from the experimental ones. 

Once the conductance of graphene is computed, the 
curves Rgr{uj) can be calculated via the Fresnel equa- 
tions using the known optical constants of Si02 and Si, 
as specified in the Appendix. The latter values are well 
known and can be further refined by the fitting of the 
reflectivity of the bare substrate Rsub{^)- 

To summarize, the conductance G{LO,Vg) and the re- 
flectivity Rgr{uj,Vg) within the presented approach de- 
pend on 9 parameters: 70, 71, 73, 74, A, U, F, T and 
VcN- We applied the non- linear Levenberg-Marquardt 
modelling routind^H to directly fit the experimental spec- 
tra. In order to speed up the iterations and improve 
the convergence, the derivatives of the reflectivity with 
respect to all adjustable parameters were calculated ex- 
plicitly using analytical formulas. One obstacle to this 
approach is that the chemical potential should be de- 
termined from the equation ([s]), which in general can be 
done only numerically. Therefore we treated /j, as a fitting 
parameter and used equation ([5]) as a rigid constraint of 
the least-square minimization. 



III. RESULTS 

Reflectivity spectra and their tight-binding 
modelling 



As it was discussed in the Section |II A the self- 
normalized reflectivity R(Vg) / R{— 20V) (from now on we 
omit the index "gr" for brevity) is the most accurately 
determined quantity, which is therefore the best suitable 
for quantitative analysis. These spectra, taken at the 
substrate temperature of 10 K, are presented in Fig|3|L 
for the whole span of gate voltages used (from -100 V 
to -1-100 V with a step of 10 V). As compared to Fig|lj;, 
the spectral resolution in this figure is diminished to 5 
meV. The spectra contain rich structure that evolves in 
a peculiar fashion as a function of the gate voltage. Such 
a complicated behavior is due to the fact that all four 
bands are involved into the electronic transitions that 
affect optical properties in the considered energy range. 

The amplitude of the structures of R{Vg)/R{-20 V) 
increases as the difference between Vg and -20 V grows. 
Therefore, it is useful to plot also the differential reflec- 
tivity spectra (Figjsb) defined as follows: 



AR 



i^,Vg) 



^ Riu;,Vg + 5V)-Riuj,Vg-5V) 
R{lu, Vg + 5Y)+ R{to, Vg~5Y) 



(6) 



Such a way of showing data emphasizes certain struc- 
tures, such as those indicated with dashed lines, and 
their gate voltage dependence. Another advantage of 
this representation is that it does not require a priori 
the knowledge of the gate voltage corresponding to the 
charge-neutral state. 

We fitted the whole set of the AR/R spectra simultane- 
ously using the tight-binding parametrization, described 



TABLE L Parameter values obtained by the least-square fit- 
ting of reflectivity spectra (fit 2). All parameters, except T 
and VcN, are given in eV. 



Parameter 


This work 


'ntr'T' -.-1 — J'iOr 
Ur 1 caicuiatioir^ 


70 


3.16 ± 0.03 


2.598 ± 0.015 


71 


0.381 ± 0.003 


0.34 ± 0.02 


73 


0.38 ± 0.06 


0.32 ± 0.02 


74 


0.14 ± 0.03 


0.177 ± 0.025 


A 


0.022 ± 0.003 


0.024 ± 0.01 


r 


0.018 ± 0.003 




T 


120 ± 15 K 




VcN 


-22 ± 1 V 





in Section [II C[ We assumed that the 70, 71, 73, 74, A, 
F, T and Vcn do not depend on Vg. Since one of our 
main goals is to detect and measure the bandgap, we 
performed fits in two different ways. First, we set the 
parameter U to zero at all gate voltages (fit 1) so that 
the difference between the spectra is only due to a vari- 
ation of the chemical potential. In the second run (fit 
2), the bandgap was allowed to vary as a function of Vg 
in such a way that U at each value of the gate voltage 
was treated as an independent parameter. In each case, 
we tried different sets of initial parameters (within the 
scope of physically reasonable values) and checked that 
the fitting routine converges to the same result. The pa- 
rameter confidence limits were estimated based on the 
correlation analysid^ and by repeating the process after 
varying data points within their error bars. 

The model curves corresponding to the fits 1 and 2 
are shown in Figures |3|l and [sJd (solid green and dashed 
red lines respectively)?^ Both fits show almost the same 
match to the data outside the region around 0.4 eV. How- 
ever, within this region the fit 1, which does not involve 
the bandgap, is qualitatively worse. It fails to reproduce 
some strong structures, in particular, the ones marked 
with the circle. As discussed in Refs ll8l20l28l this is ex- 
actly the region, where the bandgap is expected to affect 
the spectra. At the same time, the quality of the fit 2 
is remarkably good. Thus our data unequivocally show 
the presence of the bandgap. There are still some mis- 
matches that we shall address separately. 

The parameters of the fit 2, apart from U, which de- 
pends on the gate voltage, are given in the Table |T] 
One can see that the SWMcC parameters can be de- 
termined from the infrared spectra. Except 73, these 
param eters were already determined in previous infrared 
studie^i^E^E^Ilo]^ |-jy monitoring the gate voltage depen- 
dence of easily recognizable spectral features, such as the 
maximum of the 71 peak. Using the least-square fitting 
method, the parameter 73 can now also be estimated. 
This term results in a deformation of the 71 peak, but 
its effect on the spectra is more complicated than just a 
broadening produced by F (that we find to be about 18 
meV). One should mention, that in the fit 1 (where the 
bandgap was not included) this parameter relaxed to an 
artificially large value (about 0.5 eV), in order to mimic 
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FIG. 3: (color on-line) (a) Self-normalized and (b) difference reflectivity spectra as a function of Vg. Dots are the experimental 
data points, solid and dashed lines represent model fits 1 and 2 (without and with a bandgap) respectively, (c) The real part of 
the optical conductance. Curves - the result of the reflectivity flts (the colors and the line types correspond to panels a and b), 
symbols - the refined conductance obtained using a KK-consistent inversion as described in the text. In all panels the curves 
are vertically displaced for clarity. 



somehow the bandgap induced smearing of the 71 —peak. 

For comparison, we also reproduce the SWMcC param- 
eters obtained from a density functional theory (DFT) 
calculation on graphite by Charlier et (a more com- 
plete overview is given, for example in RefjlSp. The 
agreement is good, except for 70, for which the DFT 
values are somewhat lower. Nevertheless, 70 deduced 
from various experiments on graphitd^ is in an excellent 
agreement with our result. Based on the obtained value 
of 7o, we find the Fermi velocity vp = (3/2) 070 /fi. to be 
1.02 ±0.01 X lO'' m/s (a = 1.42 A is the nearest-neighbor 
interatomic distance). 

The gate voltage corresponding to zero doping, Vcn, 
can be determined with a very good accuracy. As it was 



mentioned in Section II A this value is slightly different 
from the one extracted from the simultaneous transport 
measurement. It is possible that in the case of large flakes 
and a non-optimal geometry of electrical contacts, it gives 



a more accurate value than the one given by the maxi- 
mum of the DC resistivity as it is not affected by the 
distribution of the measurement currents in the flake. 
Infrared spectroscopy can be used therefore as an inde- 
pendent indicator of the doping level. 

We flnd the broadening F to be about 15-20 meV. 
This value is larger than the one reported in RefTSI 
(w 0.0271 = 8 meV), which is perhaps due to a higher 
concentration of charging impurities (the ones that shift 
the charge neutral point from the zero bias). At the 
same time is considerably smaller than the broadening of 
about 60 meV found in Ref|22]on double gated graphene, 
which is probably related to extra scattering and/or in- 
homogeneity introduced by the top gate. 

At flrst surprisingly, the deduced temperature of 
graphene T is of the order of 100 K, even though the 
substrate was kept at 10 K. Although the true graphene 
temperature may indeed be somewhat higher due to a 



7 



weak thermal contact between the warped flake and the 
substrate, another plausible explanation is that this is an 
indication of the spatial inhomogeneity of the chemical 
potential. It is easy to see that if we neglect the change 
of the bands as a function of /i then a smearing of the 
chemical potential /i ^ ± J/i has almost the same ef- 
fect in the Kubo formula as increasing the temperature 
(ksTeff ~ (5/z). Thus we get an upper limit of about 
10 meV to the inhomogeneity of the chemical potential. 
It is worth emphasizing that the thermal and the scat- 
tering induced broadening given by T and T respectively 
are clearly distinguishable by the fitting routine. 
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FIG. 4: (color on-line) The extracted values of \ U\ (top panel) 
and the chemical potential (bottom panel) as functions of the 
gate voltage and doping. Dash-dotted line - the calculated 
"unscreened" value of |(7|, solid line - an ab initio DFT calcu- 
lation from Re£41j which takes screening effects into account. 

Figure [4|i shows the extracted bandgap as a function of 
the gate voltage and doping. Here we take \U\ as a mea- 
sure of the bandgap, since it practically coincides with 
Ag in the considered doping range. One can see that at 
small gate voltages the bandgap goes to zero within error 
bars and it grows almost linearly for both electron and 
hole doping, reaching 70-80 meV at the maximum ap- 
plied gate voltages. Interestingly, the minimum is closer 
to zero gate voltage (where n « 2 x 10^^ cm^^) than 
to Vg = VcN {n = 0). As discussed by Castro et al.'', 
such a shift can be understood by considering the dopant 
molecules adsorbed by the surface acting as an effective 
top gate electrode. However, in this case, the zero-gap 
point is expected to be at Vg = —Vcn- Seeing zero gap 
at zero gate voltage would be expected if the dopants 
are intercalated between the carbon layers, so that they 



do not introduce an interlayer electrostatic asymmetry, 
however we do not have any independent experimental 
verification of this happening. If the dopants were below 
the flake, in this picture one expects the gap to vanish at 
the charge neutral point. 



-100 V 



-20 V 



V 



1-100 V 








Momentum 



FIG. 5: (color on-line) Electronic bands (calculated using pa- 
rameters extracted from optical spectra) and chemical poten- 
tial at selected gate voltages. 

From Figure |4Jd one can see that the chemical poten- 
tial shows a monotonic, slightly sublinear increase with 
doping. In figure [5] the calculated band structures, cor- 
responding to the gate voltages -100, -20, and -1-100 V, 
are presented, together with the position of the chemical 
potential. 



B. Optical conductance 

FigjSj; shows the real part of the optical conductance 
G{uj) calculated using the obtained tight-binding param- 
eters as described above. As in FigjS^, green solid and 
red dashed lines correspond to the fits 1 and 2 respec- 
tively. The conductance is normalized by 2Go, which is 
the theoretical asymp totic value for bilayer graphene at 
high frequencied^^'^ir One can see that the model curve 
corresponding to the fit 1 shows a very broad 71 at pos- 
itive gate voltages because the parameter 73 in this fit 
was unrealistically large, as discussed above. 

Since some deviations between the model fits and 
the experiment are present, we refined optical conduc- 
tance using the Kramers-Kronig constrained variational 
methocP2. Within this approach, we represented the con- 
ductance as a sum of the two terms: G{oj) — Gmod{^) + 
Gvar{'^)- The model term is calculated using equations 
([2]), ^ and ([4|. Gvar{i^) is a variational Kramers-Kronig 
constrained correction, needed to reproduce all remain- 
ing fine details of the experimental reflectivity spectra."^. 
At the refinement stage, Gmod{^) was fixed and G^ari'^) 
was adjusted in order to get the perfect match to the 
reflectivity spectra. The refined conductivity is shown 
in Fig. [Sj: with symbols. Since we based our analysis 
on the relative reflectivity spectra, this procedure gives 
most accurately the relative changes of G{ljS) as a func- 
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tion of Vg and w (the accuracy is better than 0.1 Go), 
while the error bars for the absolute level of G{uj) can 
be somewhat larger. This explains slightly negative val- 
ues of Re G{u!) at low frequencies at high gate voltages. 
Weak structures below 0.15 eV are artefacts coming from 
the optical phonons in Si02, which are not fully cancelled 
in the fitting procedure, probably due to a weak depen- 
dence of these phonons on the electric field, which is not 
included in our model. 



C. Temperature dependence 

In addition to tuning the spectra by the gate volt- 
age, varying the temperature provides another impor- 
tant piece of information. As an example, figure[6^ shows 
Ai?(w, —A5V) I R taken at the substrate temperature Tg^b 
of 10 K, 150 K and 300 K. The spectra clearly change 
with cooling down. In general, the structures are getting 
sharper at low temperatures. However, the sharpening 
is far from being simply a uniform broadening, such as 
due the electronic scattering (parameter T). In particu- 
lar, the peak at 0.35 eV does not change, while the dip at 
0.45 eV and a structure at 0.2 both show a pronounced 
temperature variation. 
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FIG. 6: (color on-line) (a) The difference reflectivity at Vg = 
—45 V at three temperatures of the substrate: 300 K, 150 K 
and 10 K. (b) Model curves of the same quantity, calculated 
by varying the effective temperature of graphene and keeping 
other parameters unchanged (shown in Table |l|. 

Panel (b) shows calculations of the same quantity at 
the following set of graphene temperatures: 20 K, 120 K 
200 K and 300 K. In all cases, the same parameters, ex- 
cept T, were used (Table |l]). The calculated temperature 
dependence reproduces very well the experimental one if 
one assumes that the effective temperature of graphene 
is higher than Tgub- As we discussed above, this temper- 
ature mismatch may be in part due to the spatial broad- 
ening of the chemical potential. One can see that at T = 
20 K (dashed line), the spectral structures are expected 
to be much sharper than at 120 K. Similar spectra com- 
parisons at other gate voltages (not shown) provide the 
same results. 

The observation that the effect of temperature on the 



spectra is highly frequency selective is explained by the 
fact that only electronic transitions, for which either ini- 
tial or the final state are close to the Fermi level, are 
affected by the temperature. As one can anticipate from 
the good match between spectra in panels (a) and (b) 
of Figj6] the application of the same least-square fitting 
procedure at higher temperatures provides model param- 
eters, which are essentially the same as the ones pre- 
sented. Therefore we focus largely on the low temper- 
ature results. Nevertheless, the ability to quantitatively 
predict spectra at high temperatures based on the results 
obtained at low temperature corroborates the consistency 
of the used model approach. 



IV. DISCUSSION 

A. Seeing bandgap optically: zero versus finite 
doping 



ReG 

~2g7 



1 1 1 1 1 1 1 1 r- 

u=0.15eV y-peak 

' 0.5 P 




0.4 0.6 
ho)[eV] 

FIG. 7: (color on-line) Model demonstration of the effect of 
the bandgap on optical spectra in the case of finite doping and 
the case of zero doping. Top panel: the chemical potential 
is larger than the bandgap; bottom panel: p is zero. The 
inset shows the corresponding bands (energy, in eV, versus 
momentum near the K point). A simplified set of SWMcC 
parameters is used: 70 = 3 eV, 71 = 0.4 eV, 73 = 74 = A = 0. 
The broadening F is 0.01 eV. 

Optical spectroscopy is routinely used for the bandgap 
measurements in usual semiconductors due to the 
fact that only the photons, which energy exceeds the 
bandgap, are absorbed by electron-hole excitations. 
However, bilayer graphene is special in a sense that the 
bandgap is not intrinsic but is induced by the gate volt- 
age. Therefore, if only one gate is attached to the fiake. 
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as in the present case, then the gap opening is inevitably 
accompanied by doping. Making two gate electrodes (on 
top and on bottom) would induce the bandgap without 
doping. Figure [7] demonstrates, using a simplified set of 
SWMcC parameters, the effect of the bandgap in both 
cases. The bottom panel corresponds to the undoped 
case {fi = 0). One can see that an infrared absorption 
threshold appears at hio ~ Ag, due to the transitions be- 
tween the bands 2 and 3 across the bandgap (the bands 
are numbered in the inset). Such a structure was ob- 
served indeed in a recent paper by Zhang et al^, who 
used a double-gate bilayer graphene device. The second 
notable effect is a shift by Ag/2 of the second threshold 
at 71 corresponding to the transitions between bands 1 
and 3 and also to ones between bands 2 and 4. 

The top panel describes the case of a finite chemical 
potential (electron doping) . There is a striking difference 
between two cases. Now the opening of the gap does not 
produce an absorption threshold, since the transitions 
across the gap are blocked by the Pauli principle. How- 
ever, the gap affects the lincshape of the 71 peak that 
originates from a combination of interband transitions 
3^4 and 2 — > 4. Most notably, a satellite peak at 
about 71 + Ag/2 shows up. In addition, a shoulder at 
about 71 — Ag/2 appears. The satellite and the shoulder 
stem from the transitions 2 — > 4 and 3 — > 4 respectively 
close to the K point, where the are separated by the 
bandgap energy. Thus, in the doped case the only way 
to measure the bandgap is to analyze the shape of the 71 
peak. 

In reality, the position and the shape of the peak are 
also affected by the parameters 73, 74 and A, not in- 
cluded in the above demonstration, and further broad- 
ened by electronic scattering. Therefore, when the gap 
is small, its extraction from the optical spectra requires 
direct fitting of the data using a complete set of SWMcC 
parameters. When the gap is large, the identification of 
the gap becomes easier as the satellite to the main peak 
is more pronounced. We note that in actual data one can 
clearly recognize a satellite to the 71 peak for > 80 V 
(shown by arrows in Figj3|:). 

Although the fit reveals the presence of the bandgap 
also at negative gate voltages (Figj4^), the conductance 
spectra do not show a clear satellite at this doping side. 
Such a difference is in part due to the electron-hole asym- 
metry, which results in a stronger broadening of the 71 
peak at the hole doping, and in part due to the shift of 
the charge neutrality point. However, in section [IV C| we 
shall demonstrate direct signatures of the bandgap for 
both polarities of the gate voltage. 



B. Bandgap: the role of self-screening 



As discussed in section II C the bandgap is determined 
by the parameter U, which is defined as the difference 
of the electrostatic potential on the two l ayers. As it 
was extensively discussed in the literatur d^ l ^ l ^^ l '^^ l '^^ l ^^l . 



the self-screening of the external field plays a crucial 
role in the determination of U. Our data fully agree 
with this. One can see (Figure |4| that the experimen- 
tal value of the bandgap is more than two times smaller 
than the "unscreened value" (dashed-dottcd line), given 
by the external field multiplied by the interlayer distance 
(3.35 A)3ii (we assumed that the bandgap vanishes at 
Vg = 0, based on the experimental results, which means 
that charging impurities do not introduce any imbalance 
of the interlayer potential). The same observation was 
made in Refill] 

A proper microscopic calculation of the bandgap must 
be done self-consistently, since the screening depends 
on the bandgap and vice versa. Such a complicated 
problem was treate d on the Hartree level based on the 
tight-bi nding m ode F*''*^^ ^ as well as using ab-initio 
methodaSHMl] These calculations provide the doping 
dependent bandgap, which is much closeJ^to the present 
experiment than the "unscreened" model. As an exam- 
ple, we present on Figure |4]the ab-initio DFT calculation 
of Cava et (solid line), which shows a good agree- 
ment with the experimental data. 



C. "Photon energy - gate voltage" mapping of the 
interband transitions 

We saw that the extraction of the optical conduc- 
tance from the measured spectra is rather involved in 
the present case, where the measured spectra depend on 
both real and imaginary parts of G{uj) (as detailed in 
the Appendix). Now we propose a simple way to visu- 
alize electronic transitions based on the raw reflectivity 
data, which most clearly demonstrates the electron-hole 
asymmetry, the opening of the bandgap and other feature 
of the band structure. 

In figure |8^, the whole set of experimental spectra 
AR/R is represented as a color map in the coordi- 
nates {hio\Vg). One can see a set of lines that resem- 
ble somewhat band dispersions seen in ARPES. First we 
note the two "<"-like structures, shifted with respect to 
each other along the photon energy axis. They corre- 
spond to the onset-like features in the optical conduc- 
tance, marked in Fig [7] as A and B and related to the 
interband transitions 2 — > 3 and 1 — > 3 respectively 
(2^3 and 2 ^ 4 in the case of hole doping). In- 
deed, one can see that they match closely the expected 
threshold energies lj A{Vg) = 2|/x(l/g)| (dashed line) and 
^B{Vg) = 71 -|-2|/i(Vg)| (dotted line). Although several 
experimental p apers presented infrared spectra of gated 
bilayer graphen d^^ l ^^ l ^^ l ^^ l ^^ l ^^ l, the second threshold was 
reported only in Ref|3D] Here we reaffirm, based on a new 
set of data, the existence of the second threshold, which 
is essential for the overall consistency of the tight-binding 
approach. 

In figure [8^ the presence of electron-hole asymmetry 
is quite obvious, since in the case of perfectly symmet- 
ric bands with respect to the Dirac point the AR/R 
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FIG. 8: (color on-line) Color maps of [AR/ R)(ijj,Vg). (a) - experiment, (b) - model curves corresponding to Fit 1 (zero 
bandgap), (c) - model curves corresponding to Fit 2 (including bandgap), (d) - model curves in the absence of the bandgap 
and electron-hole asymmetry. The color scheme is the same in all graphs. The dashed and dash-dotted lines in panel (a) are 
the dependencies of lua = 2i/ij and lob ~ + 'i\^J^\ respectively. 



spectra should be precisely antisymmetric with respect 
to Vg ~ VcN, as exemplified in the hypothetical graph 
of Fig. [Sji. Within the SWMcC model, the asymmetry 
between electron and holes bands is due to the hopping 
term 74 and the on-site energy difference A. In Refs. 
I18ll9l thev were deduced from the doping dependence of 
the position of the maximum of the 71-peak. Since the 
maximum location is affected not only by 71, 74 and A 
but also by 73, r and, most importantly, by U , we choose 
to determine all parameters, including 74 and A, by fit- 
ting of the whole set of spectra. 

This way of presenting spectra also allows us to see 
distinct features related to the opening of the bandgap. 
These features appear to be quite different on the two 
doping sides due to the electron-hole asymmetry. In 
panels (b) and (c) of the same figure, we show the fits 
of AR/R without and with the bandgap respectively 
(namely, fits 1 and 2 described in Section III A I . On the 
electron side {Vg > Vcn) the 'ridge' indicated as '1' finds 
absolutely no counterpart in the fit 1, but is mimicked 



by in the fit 2. On the hole side, the ridge marked as '2' 
clearly disperses towards low frequencies as the absolute 
value of Vg is increasing. This trend is well captured by 
the fit 2, while in the fit 1 this ridge is precisely vertical. 



D. Deficiencies of tiie tiglit binding description 
with constant scattering 



Although the overall agreement between the panels (a) 
and (c) of Fig. [Sjis very good, a closer inspection reveals 
some deficiencies of the fit 2. For the electron doping, at 
gate voltages between 50 and 80 V, the ridge indicated 
by 1 is quite narrow in the experiment but is broad and 
barely recognizable in the fit 2 (as indicated by 3 in the 
panel (c). For the hole doping, the fit 2 contains a weak 
extra ridge (marked as '4') which is not clearly present in 
the experiment. In Figj9]we concentrate on these doping 
levels, taking the gate voltages Vg = -95 V and -1-65 V 
as examples. Here we improved the match even further 
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as compared to Figj3]3 by slightly compromising the fit 
quality at other gate voltages (although all the bandgaps 
were kept the same). 
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FIG. 9: (color on-line) The experimental curves {AR/ R){uj) 
(solid lines) at Vg = -95 V (left panel) and -1-65 V (right panel), 
together with the best fitting curves (dash-dotted lines). 

The model curve at -95 V t is featured by three dis- 
tinct peaks between 0.3 and 0.45 eV. Crudely speaking, 
they are related respectively to the shoulder, the main 
maximum and the satellite to the 71 peak depicted in 
Fig|7] The value of U is related to but somewhat smaller 
than the distance between the leftmost and the rightmost 
peaks. We can see that the data also show three peaks at 
the same photon energies, but the first peak is stronger 
and the third one is much weaker than the counterpart 
on the theoretical curve (the latter one gives rise to the 
feature '4' in Fig|9];). 

At Vg =-f-65 V the peak at 0.4 eV is rather sharp in 
the experiment but has a pronounced double structure in 
the model. Although this does not question the existence 
of the bandgap as such (recall that the best fit without 
the bandgap does not show this peak at all), such a dis- 
crepancy is too significant to be ignored. 

Presently, the origin of the shown mismatches is not 
clear. This may be an indication that more hopping 
terms need to be taken into account. In this respect, it 
would be instructive to compare optical spectra directly 
to the results of ab-initio band structure calculations. 
Another possibility is that the discrepancies are caused 
by our assumption that the scattering F is the same for 
all electronic states, which can fail due to the electron- 
phonon and electron-electron interactions. Last but not 
least we have assumed in our analysis a rigid band model, 
i.e. the tight binding band parameters (apart from the 
bandgap) are assumed to be independent of doping and 
gate voltages. That this may not be strictly the case was 
experimentally shown in an ARPES study on epitaxial 
graphene^, where an increase of 71 by about 3 % was 
observed, when U changed from to 100 meV. Electron- 
correlation effects introduce a doping and energy depen- 
dent renormalization of the bare dispersion. Also the 
gate voltages infiuence the interatomic tunneling matrix 
elements, which in turn affect the tight binding param- 
eters. Studying the manifestation of these interactions 
in optical spectra will undoubtedly be one of the most 
intriguing directions in the further research of graphene. 



V. SUMMARY 

We presented a detailed analysis of infrared reflectiv- 
ity spectra of bottom gated bilayer graphene that al- 
lowed us to determine the tight-binding Slonczewski- 
Weiss-McClure parameters and the doping dependence 
of the bandgap induced by the electric field generated by 
the gate. The direct least-square fitting of the whole set 
of infrared spectra using the SWMcC Hamiltonian and 
the Kubo formula turns out to be a very efficient tech- 
nique to disentangle the complicated interplay of vari- 
ous band structure parameters in the optical spectra. It 
also provides independent information about the extrin- 
sic doping level. 

Our analysis clearly shows the presence of the bandgap, 
which depends almost linearly on the gate voltage. This 
dependence agrees with the tight-binding and ab initio 
calculations that take the screening of the external field 
by the tt bands into account. At the maximum applied 
gate voltage of 100 V the bandgap reaches about 80 meV, 
which is three times larger than fc^T at room tempera- 
ture. Even higher values of the bandgap (up to 250 meV) 
could be obtained on double-gated bilayer graphene^, 
making this material very promising for applications. 

The very fact of achieving quantitatively good fits is a 
strong indication that the tight-binding model is quite ac- 
curate for the actual band structure of bilayer graphene. 
Nevertheless, some discrepancies remain, and further in- 
vestigations will be needed to explore their origin in the 
context of electron-phonon and electron-electron interac- 
tions. 
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F. Mauri, E. McCann, L.A. Falkovsky and D.N. Basov 
for illuminating discussions. This work was supported 
by the Swiss National Science Foundation (SNSF) by 
the grant 200021-120347, through the National Center 
of Competence in Research "Materials with Novel Elec- 
tronic Properties-MaNEP" . 



APPENDIX: Relation between the reflectivity and 
the optical conductance of graphene 

The refiectivity of bare substrate, and graphene on top 
of the substrate (Fig|lji) can be calculated based on the 
optical conductance of graphene G{uj) and the known di- 
electric functions e{uj) of Si02 and Si^°. We can treat 
the silicon layer as semi-infinite, since in our case it is 
thicker than the penetration depth. The Fresnel equa- 
tions for the refiectivities can be written as follows: 



Rs 



^01^100 

1 - riori2(/'2 

^01^100^ 

1 - fwri2(jP' 



(7) 

(8) 



where indices 0,1 and 2 refer to vacuum (e = 1), Si02 
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and Si layers respectively and 



exp 



(9) 



We used the complex reflection and transmission coefli- 
cients at the interface between media i and j: 



(10) 
(11) 



The presence of graphene between layers (in our case 
it is between vacuum and Si02) modifies the interface 
coefficients in the following way: 



2^/e 



(12) 
(13) 



<J0 



where a — /he is the fine structure constant. The latter 
formulas are valid in the thin-film limit (the thickness is 



much smaller than the wavelength), which is perfectly 
applicable to graphene. 



Since the typical values AR/R (FigjSp) are rather 
small (^ 10~^), it is useful to introduce the so-called 
" sensiti vity" functions that we previously used in similar 
analysed^^EH and employ an approximate linear relation: 



AR{uj) 
R 



/3i('^) T^^ + P2{u;) (14) 



Go 



Go 



Here we obtained the sensitivity functions /3i(w) and 
f32{uj) numerically, using a linear regression of the ex- 
act formulas for the values of G/Gq ~ 1. These function, 
which are specific to the substrate used are shown in 
(Fig 10 1. One can see that the refiectivity depends on 



both the real and the imaginary parts of G(tj) in a non- 
trivial way. At high energies, the dependence is stronger 
than at low energies. Optical phonons in Si02 give rise 
to structures at ^ 0.15 eV, which affect AR/R. 
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